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Abstract. A trivalent diagram is a connected, two-colored bipartite graph (par- 
allel edges allowed but not loops) such that every black vertex is of degree 1 or 3 
and every white vertex is of degree 1 or 2, with a cyclic order imposed on every 
set of edges incident to to a same vertex. A rooted trivalent diagram is a trivalent 
diagram with a distinguished edge, its root. We shall describe and analyze an 
algorithm giving an exhaustive list of rooted trivalent diagrams of a given size 
(number of edges) , the list being non-redundant in that no two diagrams of the 
list are isomorphic. The algorithm will be shown to have optimal performance in 
that the time necessary to generate a diagram will be seen to be bounded in the 
amortized sense, the bound being independent of the size of the diagrams. That's 
what we call the CAT property. One objective of the paper is to provide a reusable 
theoretical framework for algorithms generating exhaustive lists of complex com- 
binatorial structures with attention paid to the case of unlabeled structures and 
to those generators having the CAT property. 
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0. Introduction 



Roughly speaking, a trivalent diagram is a connected graph with degree conditions 
on its vertices and cyclic orientations on the edges adjacent to each vertex It is the 
combinatorial description of an unembedded trivalent ribbon graph [2n 113] (c/. 

for a precise definition). We shall see (c/. theorem [1.1 ) 
that it can be described by a pair of permutations {a,, do) satisfying the conditions 
of involutivity al = id and triangularity = id. The notion of rooted trivalent 
diagrams is also very useful, both to our study and to the target applications; so we 
take a special care to study them in detail. 



definitions [1.1| and 1.2 



0.1. Problem Statement. We shall describe and analyze two algorithms, the first 
giving an exhaustive list of rooted trivalent diagrams of a given size [cf. definitions 



1.3 and 1.4 below) and the second giving an exhaustive list of unrooted trivalent 
diagrams (definitions 1.1 and 1.2), those lists being non-redundant in that no two 
diagrams of the lists are isomorphic. The algorithm for rooted diagrams will be 
shown to have optimal performance meaning that the time necessary to generate a 
diagram is bounded in the amortized sense. What is striking is that the bound is 
independent of the size of the diagrams being generated. One objective of the paper 
is to provide a reusable theoretical framework for algorithms generating exhaustive 
lists of complex combinatorial structures with attention paid to the case of unlabeled 
structures and to those generators having the CAT property. 



0.2. Motivations. In a recent paper j2l] we gave a complete classification of the 
subgroups of the modular group PSL2(Z) and their conjugacy classes by rooted 
trivalent diagrams and trivalent diagrams. A question one may ask is how to gener- 
ate a complete list of such trivalent diagrams. Such a question is unavoidable: for a 
classification to be fully satisfactory there should be a systematic way to enumerate 
all the particular instances of the objects being classified. Moreover, it was soon re- 
alized that there is a connection with combinatorial maps. In this paper we clarify 
that point and give as an application a way to generate exhaustive lists of triangular 
combinatorial maps. 

The other sources of motivations to generate trivalent diagrams come mainly from 
mathematical physics in connection with two-dimensionnal quantum gravity and 
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the Witten-Kontsevich model . Algebraic topology is also a source of motivation 
through triangular subdivisions of surfaces, knots, braids, links and tangles theory 
pn [2] . It is also connected to the deformation theory of quantized Hopf algebras 
[3 [H] . The problem we solve is also relevant to the study of combinatorial maps as 
explained in section [s] and to the vast galoisian program of A. Grothendieck [9] as 
explained in hundreds of papers and books such as IT5| [6] . As an application, 

we give in section |4] a way to generate a complete list describing all the sub-groups 
of a given finite index in the modular group PSL2(Z) and a way to decide conjugacy 
relations among those subgroups. We show also in section |5} as a second application, 
how to generate an exhaustive list of triangular maps satisfying various criteria. 

0.3. What is a CAT Generator. The expression "CAT generator" is an acronym 
for constant amortized time generator meaning a generator of combinatorial struc- 
tures that average spend only a constant time generating each of the structures. 
The usual idea in such a generator is that passing from one structure to the next 
requires only a few modifications to be made. Sometimes, though, it could take 
more modifications than usual and we don't usually have any upper bound on the 
amount of actual modifications that could be needed to pass from a structure to the 
next. When need for large number of modifications tends to be significantly rare 
in comparison to small ones, we can sometimes prove that an amortization effect is 
going on. Technically, one can summarize that amortization effect in saying that 
the total amount of time needed to generate n distinct structures is asymptotically 
bounded by a constant multiple of the number n of structures being generated, the 
word constant meaning that the bound is independent of the size of the structures 
being generated. 

1. Trivalent Diagrams 

Definition 1.1. A trivalent diagram is a connected, two colored bipartite graph 
(parallel edges allowed but not loops) such that every black vertex is of degree 1 
or 3 and every white vertex is of degree 1 or 2, with a cyclic order imposed on the 
edges incident to each vertex. The size of a trivalent diagram is the number of its 
edges. 

Given a trivalent diagram F, we denote by r_, F, and To the sets of its edges, 
black vertices and white vertices respectively. Given an edge a G r_, we denote 
by d,{a) G F, and do{a) G Fq the black vertex and the white vertex to which it 
is incident. Given an edge a G F_, we denote by a, (a) and cro(a) the next edge 
incident to d,{a) and do{a) respectively, in the cyclic order. According to the degree 
conditions of the definition we have o"^ = cTq = id which implies that both o", and 
(To are bijections, so they are permutations on the set of edges F_ of the diagram. 
The connectivity condition of the definition is equivalent to the transitivity of the 
permutation group generated by a, and cxo. 

Definition 1.2. A morphism (f between two trivalent diagrams F and F' is a triple 
of mappings ip,, ipo and from the three sets F,, Fq and F_ to the three sets F',, 
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Figure 1. A trivalent diagram is conveniently described by a little diagram 
like the one above, hence the name. The actual cyclic orientation of the 
vertices are conveniently rendered implicit by adopting the counter-clockwise 
orientation of the figure. 



and r'_ respectively, compatible with the three structure mappings and the two 
permutaion in that </?_ is both equivarient to the a, and ctq permutations and the 
following diagram is commutative. 



a. 



r_ 



do 



iPo 



V 



d, do 

When those three mappings are bijections the morphism is an isomorphism. 



A important fact, well known to the experts is recalled in the following theorem. 
It is used throughout the article to formulate the algorithms. 

Theorem 1.1. The set V_ and the two permutions a, and ctq altogether suffice 

to describe the isomorphism class of the digram T. Moreover, the cycles of the 
permutations cr, and (To are in natural bisection with the black and white vertice of 
r respectively. 

Proof. Given a trivalent diagram F, an isomorphic trivalent diagram V can be 
reconstructed from the set r_ and the two permutations cr, and cXo of r_. Let 
r_/(T,, the cyles of the permutation u,, be its set of black vertice and let r_/(To, 
the cycles of the permutation do, be its set of white vertice, and define its boundary 
mappings &, and d'^ to be the natural projection of the quotients. 

We now construct the isomorphism the following way. Since d, and do are equiv- 
arient to the permutations cr, and ctq respectively, they induce natural maps </?, and 
</7o completing the following commutative diagram. 

d' d' 

r: ^ — r'_ — r; 



Taking ip_ to be the identity mapping, one has a morphism from the diagram V to 
the diagram F. To show it's an isomorphism one has to show the bijectivity of the 
three mappings ip,, ipo and The mapping being the identity, is necessarily 
bijective. The bijectivity of the mappings y?, and means that two edges are in 
the same cycles of the respective permutations cr, and o"o if and only if they are 
incident to the same black and white vertex respectively, which is guaranteed by the 
definition. □ 

1.1. Rooted Trivalent Diagrams. The following notion play an important role 
in that dissertation and in the applications. 

Definition 1.3. A trivalent diagram is said to be rooted if one of its edges is dis- 
tinguished from the others as its root. 

A convenient way to describe the rooting of a diagram is to draw a little cross on 
its distinguished edge. 

Definition 1.4. A morphism Lp of rooted trivalent diagrams (F, a) and (F', a') is a 
morphism of the underlying diagrams (ignoring the roots) which ip_ component is 
further assumed to send root to root. 

1.2. Labeled versus Unlabeled Diagrams. Historically, the dichotomy between 
labeled and unlabeled structures had been greatly clarified and properly emphasized 
by the introduction by A. Joyal of combinatorial species [TT]. The subject was, and 
still is, a very prolific source of discovery from the Quebec school of combinatorics 
and from a growing community of researchers around the world. One must cite 
the wonderful book [31 S] by F. Bergeron, G. Labelle and P. Leroux, which gives a 
briliant exposition of the whole subject. 

On a given set of vertices X one can build different trivalent diagrams and rooted 
trivalent diagrams. We denote by D^{X) and D*{X) the corresponding sets of 
structures, we call X the labeling alphabet and we refer to diagrams one can build 
on that set as diagrams labeled by X. Any bijection q between two finite sets 
X and Y induces a bijection between the sets D^{X) and D^{Y) of trivalent 
diagrams labeled by X and Y , respectively. This induced bijection is the relabeling 
operation from D?j{X) to D^{Y). It is also referred as a transport of structure along 
the relabeling bijection q. The same considerations also apply to rooted trivalent 
diagrams and in fact to any labeled combinatorial structures. 

The above discution leads to the consideration of the Joyal Functors and 
of the two combinatorial species of trivalent diagrams and rooted trivalent di- 
agrams, respectively. In the formalism of Joyal, two labeled structures are said to 
be conjugate or isomorphic if they coincide modulo the relabeling operation. An 
unlabeled structure is then just a conjugacy class of labeled structures. We denote 
by D^{n) and D'{n) the sets of unlabeled trivalent diagrams, unrooted and rooted 
respectively, and by D^{n) and D*{n) the sets of trivalent diagrams, unrooted and 
rooted respectively and labeled by the set { 1, . . . ,n}. The symmetric group 

acts via relabeling on the set of structures labeled by { 1, . . . , n } and the sets D^{n) 
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We call the corresponding natural projections, 



7r„ : D3{n) Dz{n) 
TT, : Dl{n) ^ Dl{n) 
the condensation mappings of the combinatorial species and D\. 



2. Characteristic Labeling 



A characteristic labeling is the choice of a unique representative in every conjugacy 
class of structures. In other terms, a characteristic labeling can be seen as a natural 
section to the condensation mapping tt. The characteristic labelings that we like are 
those which are computable. We like them even more if there is an efficient way to 
compute them. 

Rooted trivalent diagrams have the enjoyable property of possessing many char- 
acteristic labelings that are computable by means of efficient algorithms. This sit- 
uation is to be contrasted with that of general graphs. No algorithm is known to 
decide in polynomial time whether two given graphs arc isomorphic, and having an 
efficient algorithm computing characteristic labeling of general graphs would render 
that particular problem trivial. What makes trivalent diagrams particular in that 
respect is not so much that they are trivalent but more in that their edges are cycli- 
cally oriented at their vertices. Indeed, general graphs with only trivalent vertices 
still suffer from the above problem. 

What we give now, is a succinct description of an algorithm producing a charac- 
teristic labeling of rooted trivalent diagrams F and having linear time complexity in 
the number of edges of F. The idea is the following: build a rooted planar binary 
tree T by depth-first traversal (in prefix order) of the edges of the diagram (not the 
vertices, we insist on the edges). Given a particular edge a of F, the two directions 
that are explored from it, are given by the two operations c, and do on the set 
of edges. We take care never to revisit a previously visited edge and we label the 
edges of F by numbers from 1 to n according to the order of their appearance in the 
depth-first traversal. 

2.1. Implementation. We need as global data an integer c and the following seven 



arrays: 



visited : X 



Bool 



4 

4:{1, 
So, Si: X - 
to, ^1 : { 1, 



{1, 

,n} 



,n} 



X 



X 



,n} 
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Figure 2. If one gives as input to the relabeling algorithm (algorithm 2.2) 
the rooted diagram shown on the left with an arbitrary initial labeling on the 
arbitrary alphabet X = {a,b,c, d, e, f,g,h,i}, it produces the characterstic 
relabeling shown on the right with numbers from 1 to 9 according to the 



depth-first traversal order of algorithm |ZT} Note the natural cutting between 
the edges labeled by 2 and 8 that arises from the depth-first traversal. 



Algorithm 2.1 wliicli is an auxiliary recursive program computes the transport bijec- 



tions £o and ii. Algorithm 2.2 is the main entry point of the relabeling process. It 
does the initialization job (line 2 to 4) and the actual relabeling of the input diagram 
(line 6 to 8). It takes as input a trivalent diagram labeled with the elements of the 
set X and rooted by the element x of X. The arrays Sq and Si and the element 
X E X are descriptions of the input diagram via its associated two permutations a. 



and (To (c/. theorm 1.1). The output diagram is encoded by the two arrays to and 
ti in the very same fashion. The visited array is used to remember the positions al- 
ready visited by the relabeling process. The integer c serves as a counter to label the 
vertices in the order they are encountered, io and ii are internal arrays describing 
the mutual inverse transport bijections between the input diagram and the output 
diagram. 

Algorithm 2.1: VISIT (x : X) 

1 begin 

2 if visited [x] then return 

3 visited [ x ] <— true 

4 io[x]^C 

5 £i [c] <— X 

6 C ^ C + 1 

7 Visit (sq [x]) 

8 Visit (si [x]) 

9 end 



2.2. Correctness. The idea behind that algorithm is quite simple and presents no 
difficulty excepting the actual proof of the relabeling being characteristic. There 
are two ways to do the proof ; one is conceptual by nature and the other is more 
technical. The particular description of the algorithm is itself part of that former 
argument. We shall give both arguments because preferring one or the other is 
simply a mater of taste. Let's give the conceptual argument first. 
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Algorithm 2.2: RELABEL {x : X) 



1 

2 
3 
4 

5 
6 
7 
8 



begin 

c ^ 



for i G X do 

|_ visited [i] /a/se 

Visit (x) 

for k E {1, ...,n} do 

ti[k]^£o[si[ii[k]] 



9 end 



One could have taken the input diagram to be labeled by the set { 1, ...,n} and 
then shown that the output labeled diagram remains unchanged if one conjugates 
the input labeled diagram according to any permutation of the labeling set. Such a 
proof would typically look rather technical if not difficult. Instead, one can rather 
abstract the labeling alphabet of the input diagram to be an arbitrary n elements set 
X, this requirement being the only assumption made on X. In particular, we make 
absolutely no assumption on its elements or on any structure that it may carry. 

A moment's thought may convince the reader that abstracting the input label set 
to X and making no assumption whatsoever on its elements indeed guarantees the 
required invariance. As this argument is a bit subtle and may seems a hand-waving 
argument to most people ; we now give another proof avoiding such considerations. 



Theorem 2.1. Algorithm 2.2 produces a characteristic relabeling of the connected 
rooted trivalent diagrams of size n - that is, to o.iT'd ti are invariant under any 
bijection from X onto another set X' . 

Proof. Any bijection g between two sets of input labels X and X' induces a con- 
jugacy of the two input permutations sq and si of X yielding two permutations 
Sq = g - So - and s'i = g- si ■ g~^ of X' . Now, putting ^o{x) = c and i'o{x') = c ac- 
cording to line 4 of Algorithm 2.1 with x' = g{x) and varying x yields I'q = Iq ■ g~^. 



Similarly, considering line 5 of the same algorithm, we get = g ■ ^l■ The per- 
mutations ^0' ^1 ^^"^ ^'i verify the following identities (by line 6-8 of algorithm 



2.2): 



tQ = ^Q■so■^l, ^0 = ^0 ■ 4 



ti = io-Si-ii, t[ = i[ 



'0 ^1? 
"'I ^1' 



and substituting i[, s'q, and s[ for their above values yields, a cancellation of the 

^0 = (4 ■ -{q- So - g~^) ■ (g ■ ii) = to, 

t[ = {io ■ Q'^) ■ {q ■ si ■ g-^) ■ {g ■ ii) = t^, 



thus proving the required invariance of the output. 
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□ 



3. Generating Algorithm 



Let's imagine that while exploring a particular rooted trivalent diagram using 
algorithm 2^ of section |2} we output a sequence of events describing the particular 



cycles of the permutations to and ti we encounter at each stage of the traversal. 
Those events could typically say for example: there we reach a new unforeseen 
black vertex (forward connection) and we label its adjacent edges c, c + 1, c + 2, or 
there we reach a previously visited white vertex (backward connection), or there we 
reach an unforeseen white vertex, etc... 

One can easily convince oneself that such a sequence of events, relying only on the 
execution march of the algorithm and not on the particular labeling of the input di- 
agram, is in fact characteristic to the diagram. If sufficiently detailed, that sequence 
of events can be used to unambiguously characterize rooted trivalent diagrams. The 
idea now would be to consider a rooted planar tree with leaves labeled by rooted 
trivalent diagrams and with edges labeled by events in such a way that the sequence 
of events one gets along any branch from the root to a leaf is the very sequence of 
events that unambiguously characterize the corresponding rooted trivalent diagram. 

We now obtain a usable principle of generation if we require two further proper- 
ties: exhaustivity, meaning that every isomorphism class of rooted trivalent diagram 
gets represented on a particular leaf of the tree and non-redundancy, meaning that 
every such isomorphism class gets represented just once. Assuming that we spend 
only a constant time on each node of that tree and that the number of those nodes 
is linearly bounded by the number of its leaves, this would provide a constant amor- 
tized time algorithm to generate rooted trivalent diagrams. 

To ease the memory requirements of the generator, we won't actually build the 
generation tree in memory. It will instead be realized in the calling pattern between 
the procedures of the generating program. Also, the program would be more useful 
if it generates the diagrams in permutational form instead of a sequence of events 
describing it. This means that we have to carry around a partial diagram that gets 
built while exploring the generation tree, each generating event completing that 
description and each backtrack reversing the particular changes we have made. 



3.1. Implementation. The generating algorithm uses as global data two integers 
c and n, a stack of integers and two integer arrays 

representing the rooted trivalent diagram being constructed by its black and white 
permutations cr, and cTq respectively. The integer n represents the maximum size 
of the diagram being generated while the integer c is the labeling counter used to 
attribute integer labels to the edges of rooted trivalent diagrams being generated. 
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Figure 3. Overall structure of the generating algorithm. 



The manipulation of the stack is done through the following five primitives. 

Push : Integer x Stack Stack 
Pop : Stack —>■ Integer x Stack 
StackIsEmpty : Stack Bool 
Mask, Reveal : Integer x Stack Stack 

The stack can be implemented using a doubly linked circular list represented by two 
zero-based arrays of integers. 

iV,P:{0,...,n}^{0,...,n} 

The item of index zero is just a sentinel and the stack is considered empty if the 
following relation holds, 

N[0] = P[0] = 0. 

The Mask and Reveal procedures implement removal and insertion primitives using 
a trick popularized by Knuth [12] under the name of "dancing link" . Namely, a call 
to the Mask procedure with parameter s removes the item s of the stack using the 
following two instructions, 

N[P[s]]^N[s], 
P[N[s]]^P[s], 

while a subsequent call to the Reveal procedure with parameter s restores the pre- 
vious state of the stack, before the call to the Mask procedure, using the following 
two instructions, 

N[P[s]]^s, 
P[N[s]] ^ s. 



The generating algorithm is composed of seven procedures (algorithm 3.1 to 3.7) 
and a user defined procedure called Output that serves as an outlet to the algorithm 
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and that can be used for instance to do printing jobs or to collect some statistics on 
rooted trivalent diagrams. The overall structure of the calling pattern between those 
procedure is shown in figure [3] The algorithm works by a recursive exploration of 
the structure being constructed in a way that mimics the depth first traversal of the 
labeling algorithm of section 2. 



The recursion has two base cases that are produced by the Generate procedure 
(algorithm 3.1) which is the main entry point of the algorithm. The inductive step 
of the recursion corresponds to a call to the Dispatch procedure (algorithm 3.2), 
which purpose is to successively handle each of the various cases one can encounter 
at each stage of the construction/exploration of the diagrams. This is the branching 
part of the generating algorithm in the sense that it is there that the generation tree 
forks into subtrees eventually leading to the leaves where the produced structures 
reside. A call to the Dispatch procedure results in a call to the Recurse procedure 
(algorithm 3.7) through each of the four Try procedure (algorithm 3.4 to 3.5). The 



purpose of the Recurse procedure is to call the Output procedure if the stack is 
empty meaning that the exploration/construction is finished and that we can thus 
output a finished structure or to pop an edge and call the dispatch procedure if the 
stack is not empty meaning that the structure being explored/constructed is not yet 
finished. 



3.1.1. The Generate procedure. The Generate procedure (algorithm 3.1) is the main 
entry point of the program. It is responsible for the two base cases of the induction, 
namely whether the produced structure has a univalent black vertice adjacent to its 
root edge (handled line 3 to 5 of the procedure) or a trivalent one (handled line 7 
to 13). 

In the case where the vertex is univalent, the corresponding fixed point is built 
(line 4 of the algorithm), the labeling counter is set to 2 (the label of the next 
encountered edge), and then the exploration/construction continues in the direction 
of the white vertex by calling the Dispatch procedure with parameter 1 (line 5 of 
the algorithm). 

In the case where the vertex is trivalent, the three edges around it are labeled 1, 2 
and 3 in counterclockwise direction, the corresponding cycle in the black permutation 
is built by the three instructions line 8 to 10 of the algorithm and the two edges 
1 and 2 are pushed on the stack for further exploration (line 11 and 12) while the 
edge 3 is explored in direction of its white vertex by calling the Dispatch procedure 
with parameter 3 (line 13 of the algorithm). 



3.1.2. The Dispatch procedure. This is the start of the induction step of the gener- 
ation algorithm. The hypothesis at its start is that the two arrays sq and si refiects 
the structure of a partial trivalent diagram being explored according to the depth 
first traversal of the labeling algorithm of section |2j The current edge s and the 
labeling counter c refiect the stage of the exploration. The exploration/construction 
is supposed to continue from the current edge s in the direction of its white vertex 
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Algorithm 3.1: GENERATE () 



begin 

it n > 1 then 

2 

so[l] ^ 1 
Dispatch 

if n > 3 then 

so[l] ^ 2 
So [2] ^ 3 

So [3] ^ 1 
Push (1) 
Push (2) 
Dispatch (3) 



14 end 



and c is the label attributed to the next unlabeled edge we encounter. At the begin- 
ning of the procedure, we don't know if that white vertex is univalent or bivalent, 
and if it's bivalent, we don't know which is the other edge incident to it. There are 
four possible cases: 



Case 1. The white vertex incident to the current edge s is univalent. That case is 
handled by a call to the TryClosedWhite procedure line 3 of the Dispatch procedure 
(algorithm 3.2). We can then assume for the three other cases that this vertex is 
bivalent. 



Case 2. The edge adjacent to the current edge by its bivalent white vertex hasn't 
been visited yet and the next black vertex is trivalent. This case is handled line 4 
by a call to the TryForward procedure. 

Case 3. As in the previous case, the adjacent edge hasn't been visited yet but 
here the next black vertex is univalent. This case is handled line 5 by a call to the 
TryClosedBlack procedure. 

Case 4. The edge adjacent to the current edge by its bivalent white vertex has 
already been visited and thus already has a label, which we call t. The point is 
that those edges, already visited but still waiting for further exploration on their 
white side, are precisely those that are stored in the stack. Each edge stored in the 
stack corresponds to an admissible possibility for the white neighbour of the current 
edge. The exploration of those possibilities is done by the "for" statement line 6 
to 9 of the procedure. The edge of the stack t matched with the current edge s is 
temporarily removed from the stack using the "dancing link" trick implemented by 
the Mask and Reveal primitives called line 7 and 9 of the procedure. 
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The production of new edges has to be compatible with the maximum allowed 
size of the diagrams n. That condition is checked by the two "if" statements line 4 
and 5 of the procedure. 

Important. We claim that those four cases cover all the possibilities and that they 
are mutually exclusive. 

Remcirk. The order in which the cases are handled by the Dispatch procedure only 

affects the order in which the structures are produced but not the way they are 
labeled nor does it change the set of structure that is produced. 



Algorithm 3.2: Dispatch (s : integer) 



1 local t : integer 

2 begin 

3 TryClosedWhite (s) 

4 ifc + 3<n + l then TryForward (s) 

5 ifc+l<n + l then TryClosedBlack (s) 

6 for t e Stack do 

7 Mask {t) 

8 TryBackward {s,t) 

9 Reveal {t) 

10 end 



3.1.3. The TryClosedWhite procedure. It handles the case where the current edge 
s is incident to a univalent white vertex (case 1. above), as the following picture 
shows. 



6' 

► — o 



It simply builds a fixed point on s in the white permutation (line 2 of the procedure) 
then it calls the Recurse procedure. 

Algorithm 3.3: TryClosedWhite (s : integer) 

1 begin 

2 Si [ S ] ^ S 

3 Recurse () 

4 end 



3.1.4. The TryForward procedure. Its purpose is to handle the case where the current 
edge s is incident to a bivalent white vertex followed by a trivalent black vertex (case 
2. above) as the following picture shows. 
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The edges incident to that trivalent black vertex are supposed to never have been 
encountered before and are then labeled by the value c, c + 1 and c + 2 (the case 
where the adjacent edge has already been encountered and thus has already a label 
is handled by the TryBackward procedure). Lines 2 to 6 build the corresponding 
black and white cycles in the permutation arrays. Among the three created edges, 
two need further exploration on their white side so they are put on the stack by the 
Push instructions line 7 and 8. Before calling the Recurse procedure the labeling 
counter c is increased to account for the creation of the three new edges. The sate 
of the stack and the labeling counter are both restored to their previous value by 
the instructions line 11 to 13, before the procedure exits. 

Algorithm 3.4: TryForward (s : integer) 
1 begin 



2 


So [c] ^ c+ 1 


3 


so[c+ 1] ^ c + 2 


4 


so[c + 2] ^ c 


5 


Si [s] ^ c 


6 


Si [c] ^ s 


7 


Push (c+ 1) 


8 


Push (c + 2) 


9 


c <— c + 3 


10 


Recurse () 


11 


c <— c — 3 


12 


Pop 


13 


Pop 


14 end 



3.1.5. The TryClosedBlack procedure. It handles the case where the current edge s 
is incident to a bivalent white vertex followed by a univalent black vertex (case 3. 
above), as the following picture shows. 



It builds the white 2-cycle and the black 1-cycle line 2 to 4 then the labeling counter 
c is increased line 5 to account for the creation of the new edge labeled c. It calls 
the Recurse procedures line 6 and then restores the value of the labeling counter 
before it exits. 

3.1.6. The TryBackward procedure. It handles the case where the current edge s 
is incident to a bivalent white vertex followed by an edge t that has already been 
visited (case 4. above), as the following picture shows. 



s 



c 



-o- 



S 
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Algorithm 3.5: TryClosedBlack (s : integer) 



1 begin 

2 Si [ S ] ^ C 

3 Si [c] 4— s 

4 So [c] ^ C 

5 C <— C+ 1 

6 Recurse 

7 C ^ C — 1 

8 end 



The edge t is chosen in the Dispatch procedure and is removed from the stack 
before the TryBackward procedure is called and reinstalled back after the procedure 
terminates. The procedure simply binds together the two edges s and t by building 
a 2-cycle in the white permutation and then calls the recurse procedure. 



Algorithm 3.6: TryBackward (s, t : integer) 

1 begin 

2 Si [ S ] 4— t 

3 Si [t ] ^ S 

4 Recurse () 

5 end 



3.1.7. The Recurse procedure. Its purpose is to check for the termination of the 
recursion. If the stack is empty then the recursion terminates and it calls the Output 
procedure. In that case, the two arrays sq and si describe a finished rooted trivalent 
diagram with edges labeled from 1, the root edge, to c — 1, the last attributed label. 
The size of the diagram is thus c — 1. Otherwise, when the stack is not empty, an 
edge is poped out of the stack and the recursion continues by a call to the Dispatch 
procedure line 7. 



Algorithm 3.7: Recurse ( ) 



1 local k : integer 

2 begin 

3 if StackIsEmpty ( ) then 

4 I Output () 

5 else 



6 
7 
8 

9 end 



k ^ Pop ( ) 
Dispatch (k) 
Push (k) 
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3.2. Correctness. In this section we show that the algorithm is correct and provide 
a complexity analysis of its execution time. This complexity analysis is based on 
a study of the structure of the execution tree of the algorithm and rely on the 
assumption that it is finite. We first prove that assumption. 

Lemma 3.1. The generating algorithm terminates in a finite amount of time. 

Proof. Looking at the procedures we see that each of them takes only a finite time 
to complete provided the procedure that are called also take in turn a finite time 
to complete. So the proof reduces to show that only a finite number of procedures 
is called, meaning that the execution tree of the algorithm is finite. Using Konig 
lemma on trees one has to show that all the branches of that tree are finite. We show 
that there is a non-negative integer quantity that is strictly decreasing along every 
branches. If we denote by the size of the stack then fi = 2ns + n — c+1 is such 
a quantity, where n is the maximum size of the diagrams being generated and c is 
the labeling counter of the algorithm. It is non-negative because is non-negative 
and because according to the "if" conditions of the Dispatch procedure, n — c + 1 is 
also non-negative, meaning that no label exceeding n is ever attributed to an edge. 
In the following table we summarize the changes in the value of n^, c and fi after a 

cycle Dispatch Try > Recurse Dispatch has been completed in each of the 

four cases described in section 3.1.2 {of. figure |3]). 





case 


1. 


case 2. 


case 3. 


case 4. 




ris - 


1 


ns + l 




Us -2 


c' 


c 




c + 3 


c+1 


c 




fi- 


2 




yU — 3 





In each case, /x' < /i so the quantity /x is strictly decreasing along every branch of 
the execution tree. □ 

Lemma 3.2. The rooted trivalent diagrams produced by the generating algorithm 
are labeled according to the characteristic labeling of section^ 

Proof. The proof is by induction. Assuming that the stack, the labeling counter 
and the labels of the already generated cycles agree with the corresponding state of 
the labeling procedure of section [2] at the start of a call to the Dispatch procedure 
(induction hypothesis) one can check that the execution of the algrithm through 
each of the four Try- procedures (algorithm 3.4 to 3.5) each preserves that hypoth- 
esis, that is the new cycles introduced in the permutations sq and si are labeled 
consistently with that of the labeling procedure and that the state of the stack also 
match the one found in the labeling after visiting those new cycles. Finally, one 
has to check that the two base cases produced by the Generate procedure are also 
consistent with the characteristic labeling. This is immediate and completes the 
induction. □ 

Theorem 3.3. The generating algorithm produces an exhaustive and non-redundant 
list of rooted trivalent diagrams. 

Proof. Exhaustivity comes primarily by induction from the local exhaustivity claim 



of the case analysis of section 3.1.2 The mutual exclusion of the cases ensures 
that different rooted diagrams are produced, at least differing in the way they are 
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labeled (the labeling counter is strictly increasing during the generation process of 
a structure), but since the structures are produced in their characteristic labeling, 
none of them could be isomorphic. □ 



3.3. Average Time Complexity. In this section we prove the main property of 
the algorithm, that it spends a constant amortized time generating each structure. 
One way to do the proof could be to express an estimate of the total execution 
time and the number of structure produced and show that the quotient of those 
two quantity is bounded independently of the size of the produced structure. We 
propose instead a proof of the majoration based on the following principle and a 
careful analysis of the execution tree of the generating algorithm. 

Balance Principle. In a finite tree, the number of leaves is greater than the number 
of its node having degree at least 2. 

Proof. The proof is by an easy recurrence on the number of internal node. Every fi- 
nite trees can be constructed starting from a one node tree by successive replacement 
of a leaf by an internal node having only leaves as sons. The starting tree satisfies 
the balance principle so the recurrence is initialized. Assuming by recurrence that 
a finite tree satisfies the principle, and replacing one of its leaves by a internal node 
having A; > 1 leaves as sons, one increases the number of internal nodes by 1 and 
the number of leaves by A: — 1 > 0. If A; > 2 the number of nodes having degree at 
least 2 is increased by 1, but the number of leaves is also increased by A; — 1 > 1 so 
the resulting tree still satisfies the principle. The recurrence is complete. □ 

Lemma 3.4. The total execution time of the generating algorithm is 0(C„), where 
Cn is the number of procedure called during the execution. 

Proof. The only procedure of the algorithm that contains a loop an thus can have 
an arbitrary long execution time is the Dispatch procedure. Since each iteration 
of the loop has a constant execution time and that each time the TryBackward 
procedure is called, we can transfer the expanse of the iteration to the TryBackward 
procedure and then assume the Dispatch procedure to have a bounded execution 
time. This way every procedure is assumed to have a bounded execution time so 
that the total execution time is proportional to C„. □ 

Let D„, Rn, On denote respectively the total number of time the Dispatch, Re- 
curse and Output procedures are called and let T„ denote the total number of time 
one of the four Try- procedures is called, n is the maximum size of the structures 
being produced. 

Lemma 3.5. We have < 2 0„. 

Proof. Let D'^^ denote the number of call to the Dispatch procedure that have an 
out degree at most 2. The leaves of the execution tree of the algorithm are calls to 
the Output procedure because the other procedures all have out degree at least 1, 
hence by the balance principle above D'^ < On- The Dispatch procedure has an out 

degree at least one because it calls the TryClosedWhite procedure unconditionally 
line 3, and when its outdegree is 1 then it means that the stack is empty (line 6) 
and that c < n (line 4 and 5). This means that all the n edges has been labeled and 
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that the call to the Recurse procedure subsequent to the call to the TryClosedWhite 
will result in a call to the Output procedure and no further call to the Dispatch 
procedure. Therefore we see that the Dispatch procedure can have an out degree 
of 1 but that can only happend once at the end of each branch. We thus have 
Dn<D'^ + On and then L>„ < 2 0„. □ 

Theorem 3.6 (CAT property). The average time spent by the generating algorithm 
to produce each structure is bounded independently of their size. 

Proof. Let C„ denote the total number of procedure called during the execution 

of the algorithm. The number of structure produced by the algorithm is equal to 
On- Since the total execution time of the algorithm is 0(C„) the average time spent 
producing each structure is 0(C„/On)- We have to show that the quotient Cn/On 
is bounded independently of n. We clearly have C„ = 1 + Dn + Rn + Tn + On, 
the +1 accounts for the first call to the Generate procedure. Since the Recurse 
procedure is only called by one of the four Try- procedures and that each one calls 
it exactly once we have T„ = Rn- Since the Recurse procedure is called twice by 
the Generate procedure and that the Recurse procedure calls one of the Recurse or 
Output procedures, we have Rn + 2 — Dn + On- Therefore we have, 

Cn^l + Dn + Rn + Tn + 0, 
^l + Dn + 2Rn + 0n 
^3Dn + SOn-S 
<90n 

and then Cn/On < 9. The bound on the quotient is independent of n as announced. 
□ 



Tn — Rn: 
as Rn + 2^ Dn + On, 

as D„ < 2 On, 



4. First Application: Modular Group and 
Unrooted Trivalent Diagrams 



We recall that the modular group PSL2(Z) is the following group of 2 by 2 integer 
matrices with unit determinant: 



PSL2(Z) = { ± J) e M2{Z)/ ± Id 



ad — be — 1 



There are many possible finite presentations for this group and we shall stick to the 
following: 

PSL2(Z) = {A,B\A'^ = B-^ = 1) 



with A and B being the following two matrices: 



-1 

1 



and B^± 



1 
-1 1 



for it renders explicit the following isomorphism: 

PSL2(Z) ~Z/2Z*Z/3Z. 
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4.1. Displacements Groups. The modular group acts naturally on the set of 
edges of any trivalent diagram. This action is generated by the two elementary 
moves, 

a ■ A = ao{a) and a-B = a,{a). 

The elementary move A acts by exchanging positions of the two adjacent edges of any 
bivalent white vertex and by fixing the only adjacent edge of any univalent white 
vertex. Similarly, the elementary move B acts by cyclically permuting the three 
edges incident to any trivalent black vertex and by fixing the only edge incident to 
any univalent black vertex. 



A 



A 




B 



Figure 4. Here is a picture of the result of the action of the two elementary 
nnoves A and B on the various sorts of edges. 



Given any trivalent diagram F, the two elementary moves just described generate 
a group $r called the displacement group of F. It is easily verified that it is the quo- 
tient group of PSL2(Z) by the kernel of the group action p : PSL2(Z) — > &r_- The 
modular group has therefore a universal status with respect to that construction, it 
can be considered as the universal group of displacements for the species of trivalent 
diagrams. If one restricts one's attention to finite trivalent diagrams, the profinite 
completion of PSL2(Z) would be a more appropriate candidate for that purpose. 



4.2. Unrooted Planar Binary Trees. One can associate to any (unrooted) pla- 
nar binary tree © a connected and acyclic trivalent diagram F, called its enriched 
barycentric subdivision F = 0''*^, by putting an extra white vertex in the middle of 
every edges of 0. The set of directed edge of and that of undirected edges of F 
are in bijection in two natural ways. 

There is another famous presentation of the modular group. It is given by two 
generators S and T and two relations as follows, 

PSL2(Z) ^{S,T\S^^ (ST)^ = 1 ). 

with S and T being the following two matrices: 

The conversion between the two presentations is done through the application of 
the following rules: 

A^S, S^A, 
B {STf, T AB-^. 
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ao 



a4 



ai3 



as 




an 



ag 



Figure 5. We see in this example the result of iterating the elementary 
move T on the edges of a binary tree. The edges are labeled by a^ where 
ttk+i = ak ■ T. This can be used to implement depth-first traversals in a 
purely iterative way. 



Here are two basic criteria relating connectedness and acyclicity of finite trivalent 
diagrams to the transitivity of the action of some displacement groups: 

1) A finite trivalent diagram F is connected if and only if its displacement group 
$r acts transitively on its set of edges. 

2) If it is a tree, then the subgroup \l/r of its displacements generated by the 
elementary move T acts transitively on its set of edges. 

There is a natural bijection between trivalent diagrams having no univalent white 
vertex and those having no univalent black vertex. It works by removing every 
univalent black vertex and the adjacent edges in one direction and by growing every 
univalent white vertex with a new edge and a new univalent black vertex in the 
other direction. This bijection is compatible with connectedness and acyclicity and 
thus restricts from the class of trivalent diagrams to that of planar binary trees, and 
we recover the classical bijection between complete and incomplete planar binary 
trees. 



4.3. Classification Principle. To any connected rooted trivalent diagram, one can 
moreover associate the subgroup of the PSL2(Z) consisting of elements that fix the 
distinguished edge of the diagram. We proved in |21] that this correspondence is 
one to one and we presented an enumeration in the form of a generating series that 



agrees perfectly with the number of structures generated by algorithm 3.1 



Now, if one changes the distinguished edge of a rooted trivalent diagram, the cor- 
responding subgroups get conjugated. We have moreover proved that two subgroups 
in the modular group are conjugated if and only if the associated rooted trivalent 
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Table 1. Trivalent diagrams of size up to five, up to isomorphism. 




Table 2. Trivalent diagrams of size six, up to isomorphism. 




Table 3. Trivalent diagrams of size seven, up to isomorphism. 



diagrams only differ by the position of their distinguished edges. It follows that the 
unrooted trivalent diagrams correspond in a one-to-one fashion to the conjuguacy 
classes of subgroups of the modular group. 
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Table 4. Trivalent diagrams of size eight, up to isomorphism. 



We gave in [21] the exhaustive hst of trivalent diagrams of size up to nine. That 
list was computed by hand in a non-systematic fashion. One intent of Algorithm 



3.1| of section |3.1| is to permit a retroactive validation of both those tables and the 
associated generating series that we reproduce here in tables |6] and [7]; they are parts 
of the online encyclopedia of integer sequences [20] under the references (A005133) 
and (A121350). 



We generated unrooted trivalent diagrams using a method already described in 
[27] : by imposing a linear order on the set of rooted trivalent diagrams with a given 
number of edges, generating them all and accepting only those that are maximal, 
according to the linear order, among all those that differ only in the position of their 
root. Table [Tj to [5] show all the unrooted trivalent diagrams with up to 9 edges. 



5. Second Application: Triangular Maps 



By an (oriented) triangular map we mean a finite polyhedral structure composed 
of vertex, directed edges, and oriented triangular faces with an incidence relation 
between them. Suggestively enough, a directed edge, also called an arc, is bordered 
by an ordered pair of vertices, which we call its origin and its destination respec- 
tively, such that triangular faces are each bordered by a cycle of three arcs whose 
destination coincides with the origin of the following arc in cyclic order. 

The following definition is useful in grasping the incidence relations of a combi- 
natorial map but insufficient because it lacks some traversal information such as the 
cyclic orientation of the faces. 

Definition 5.1. A combinatorial pre-map A is given by three sets Aq Ai and A2 
and five mappings s,t : Ai ^ Aq, £, r : Ai — A2 and : Ai — Ai satisfying the 

22 













Table 5. Trivalent diagrams of size nine, up to isomorphism. 

following conditions for all elements a of Ai, 

s{a-^) = t{a) e{a-^) = r(a) = a 

t(a-^) = s(a) r(a-^) = £(a) a"^ ^ a 

the four mappings s, t, £ and r are further assumed to be surjective. 



The elements of the three sets Aq, Ai, and A2 are the vertices, the arcs (directed 
edges), and faces of the combinatorial map respectively. The two mappings s and 

t map any arc a to its origin s{a) and destination t{a). The two mappings i and 
r map any arc a to its left-hand face £{a) and right-hand face r{a). Finally, the 
mapping .^^ maps any arc a to its inverse a^^ obtained by reversing its direction. 
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Z)«(t) = t + f + At^ + 8t^ + 5t^ + 22f + A2f + A0t^ + 120t^ + 265t^^ + 286t^^ 
+ 764 1^^ + 1729 1^^ + 2198 + 5168 + 12144 1^^ + 17034 + 37702 
+ 88958 + 136584^2° + 288270 1^^ + 682572 1^^ + 1118996 P 
+ 2306464 + 5428800 1^^ + 9409517 + 19103988 + 44701696 
+ 80904113t2^ + 163344502 f'^o + 379249288t^^ + 711598944t32 
+ 1434840718 f'^^ + 3308997062 1^^ + 6391673638 + 12921383032 
+ 29611074174F + 58602591708 f'^^ + 119001063028P 
+ 271331133136 + 547872065136 + 1119204224666 
+ 2541384297716 + 5219606253184 + 10733985041978 
+ 24300914061436 + 50635071045768 + 104875736986272 
+ 236934212877684 1^'' + 499877970985660 + o(t^^) 

Table 6. Order fifty development of the generating series Dl{t) giving as 
the coefficient of the number of connected rooted trivalent diagrams with 
n edges (A005133) which is also the number of index n subgroups in the 
modular group PSL2(Z). 



Dsit) = t + t"^ + 2t^ + 2t'^ + + 8f + 6f + 7t^ + Ut^ + 2^^^^ + 26t^^ 
+ 80 + 133 1^^ + 170 + 348 1^^ + 765 + 1002 + 2176 
+ 4682 + 6931 t^o + 13740 + 31085 + 48652 t^s + 96682 
+ 217152 + 362779 + 707590 + 1597130 1^® + 2789797t^^ 
+ 5449439 + 12233848 1^^ + 22245655 + 43480188 
+ 97330468 + 182619250 1^^ + 358968639 f'^^ + 800299302 f'^^ 
+ 1542254973 f'^^ + 3051310056 + 6783358130 + 13362733296 
+ 26648120027 f ^2 + 59101960412 1^''^ + 118628268978 
+ 238533003938 + 528281671324 1^^ + 1077341937144 
+ 2184915316390 + 4835392099548 + 9997568771074 + o{t^°) 

Table 7. Order fifty development of the generating series D3{t) giving as 
the coefficient of t" the number of connected unrooted trivalent diagrams 
with n edges (A121350) which is also the number of conjugacy classes of 
index n subgroup in the modular group PSL2(Z). 



Definition 5.2. A morphism (p between two combinatorial pre-maps A and A' is 
a triple of mappings (/jq : Aq ^ Aq, : Ai ^ A'^ and ip2 : A2 Ag compati- 
ble with the five structure mappings in the sense that the following diagrams are 
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commutatives. 






When the three mappings are bijections the morphisme is an isomorphism. 



5.1. Cyclic Orientation. In a given combinatorial pre-map A, the inner border 
of a face / is the set i~^{f) = { a G Ai | i{a) = / } of arcs having / as their left- 
hand face. A combinatorial pre-map is said to be strictly triangular if each of its 
faces has exactly three arcs in its inner border. The following definition describes the 
traversal information lacking to a triangular combinatorial pre-map to fully describe 
a triangular map. 

Definition 5.3. A (strictly) triangular map A is a strictly triangular combinatorial 
pre-map together with a cyclic order imposed on the three arcs that have the same 
left-hand face: a + 1 denotes the next arc after a in this cyclic order. 

Definition 5.4. A morphism if (respectively, an isomorphism) between two trian- 
gular maps A and A' is a morphism (respectively, an isomorphism) of the underlying 
combinatorial pre-maps that preserves the cyclic order of the previous definition. 

5.2. Associated Trivalent Diagram. To any triangular map A one can associate 
in a natural way a trivalent diagram F = A™'^ which is called the incidence diagram 
of the triangular map. To any face of A one associate a black vertex of F to any edge 
of A one associate a white vertex of F then one connect black and white vertices 
of F according to the incidence of their A counterpart. The cyclic orientations 
around black vertice is the one given by the cyclic orientation of the corresponding 
triangular faces of A. Table |8] gives three simple examples of the correspondance. 
This operation has already been described by Walsh in [26], namely: taking the 
face-vertex dual of A one gets a trivalent map and then applying the construction 
of Walsh, one gets the incidence diagram F. 

The incidence diagram of a triangular map A is the trivalent diagram, denoted 
\inc^ which sets of white vertices, black vertices and edges are the following, 

where A^ is the set of undirected edges of A, and which structure mappings do : 
^vnc _^ A*"'^, d, : A^JJ^ A*"'^ and the permutations a, and o"o of A"^^ are defined by 
the following relations, 

do{Ca) = 0^(a), C7o(Ca) = C^-l, 

d,{Ca) = bi(a), CT,(Ca) = Ca+1, 

where 7r(a) is the the undirected version of the arc a. This operation is functorial for 
it is easily extended to morphisms of triangular maps by the following process: to 
any morphism (p between two triangular maps A and A', we associate a morphism 
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denoted y?*"''^ between the corresponding incidence diagrams A*"^ and (A')™^. It is 
given by the following relations, 

Functoriality should be obvious by careful inspection. 

The .""'^ functor we get by what precedes is full and faithful, but not essentially 
surjective because trivalent diagrams we get as the incidence diagram of a trivalent 
map A have no univalent white vertex nor univalent black vertex. Let's call regular a 
trivalent diagram having no univalent vertex, the .^"'^ functor is essentially surjective 
on the full subcategory of regular trivalent diagrams, and so. 

Theorem 5.1. The functor realize an equivalence between the category of tri- 
angular maps and the full subcategory of regular trivalent diagrams. 

Proof. To prove this theorem we shall describe a reconstruction operation, which 
associate to any regular trivalent diagram F a triangular map denoted F™"^, with 
functorial property and such that for all regular trivalent diagrams F and triangular 
maps A, one have two natural reciprocity isomorphisms as follows, 

(J^mapync ^ p ^.^^ ^p^inc-^map ^ ^ 
nat. nat. 

We shall first introduce some notations. Let's call \&r the subgroup of $r generated 



by the elementary move T (c/. section 4.2) and lets call tt : F_ ^ F_/\E'r the natural 



projection. The mapping induced between F_/^r and F'_/\I^r by an equivariant 
mapping (y9 : F_ — > F'_ will be denoted y^^. The sets of vertices, edges and faces of 
the reconstructed map are the following, 

-pmap r ^ -i -pmap r -i -pmap r r -i 

i — i^x |xer_/*r J- 1 — 1 2a |aer_ i 2 — \)y Sy&. 

i he nve structure mappmgs s,t : L —>■ L ^ , r,i : L ^ 2 ■ ■ ^ 1 ^ 

F™**^ and the group action + : p™"^ x Z ^ -pmap ^^^^ reconstructed map are given 
by the following equations, 

S{ta) = ^n{a) = fo.(a) <^a^ = ^ao{a) 

t{ta) = f7r(<7o(a)) r{ta) = fa.(a-i) ta + 1 = ta,{a) 

The construction then extends to morphisms in the sense that any morphism cp 
between to regular trivalent diagrams F and F' induces a morphism ip"^"-P between the 
two reconstructed maps F™"^ and (^r')'""^ which three components are the following. 

The fonctoriality of the reconstruction operation and the two reciprocity isomor- 
phisms should be clear by careful inspection. □ 

Remark. The content of the above theorem is nothing but a specific notion of 
Poincare duality. 



5.3. Exhaustive Generation of Triangular Maps. To adapt the generator algo- 
rithm to produce only regular rooted trivalent diagrams and thus rooted triangular 
maps, it suffice to remove the call to algorithms 3^ and 3^ from line 6 and 7 of the 
Dispatch procedure (algorithm 3.2) and the lines 2 to 5 from algorithm 3.1 Those 
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Table 8. The three triangular maps with two faces together with their 
incidence diagrams. The triangular map in the middle row corresponds to 
the only way to build a torus using two triangles, and the two others, top and 
bottom rows, correspond to the two different ways to build a sphere using 
two triangles. 



Table 9. The eleven triangular maps with four faces. According to table 



15 six of them are spheres while five of them are toruses. This is easily 



checked using Euler formula. 



two removals preserve the CAT property we thus get this way a constant amortized 
time generator for regular rooted triangular maps, as announced. 



The tables [8] to 11 shows exhaustive lists of unrooted triangular maps produced 



from the output of the generator algorithm. Those tables arise from a two steeps 
process. First, generate the full list of rooted regular diagrams of a given size, then, 
remove the duplicated diagrams from the list one obtain by forgetting the root. An 
easy and efficient way to do that is described in , one has to put a linear order on 
the set of rooted trivalent diagrams of a given size, then one has to remove from the 
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Table 10. The eighty one triangular maps with six faces (first 
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Table 11. The eighty one triangular maps with six faces (last part, contin- 
ued from table 10). 



list the trivalent diagrams that are not minimal in their conjugacy class (two rooted 
diagrams being conjugated by definition, when they differ only by the position of 
their root). 



The interpretation of the drawings of tables |8] to 11 deserves some explanation. 
For that purpose we adopt a geographical terminology. Ignoring for a moment the 
surrounding circle and the dashed lines of the drawings and the triangular regions 
are called the countries of the maps and the plain lines are the boundaries of their 
adjacent countries. One can distinguish the boundaries that are bordering two 
distinct countries (the inner boundaries) from those that are bordering a single 
country (the outer boundaries). The roads of the maps are symbolized by dashed 
lines connecting, in a two-by-two fashion, the outer boundaries of the maps. 



To produce them, we have considered the planar rooted binary tree of the depth- 



first traversal of algorithm As already noted, this algorithm provides natural 
cuttings for the associated trivalent diagrams. Those cuttings arise as what we 
previously called backward connections. In contrast, the edges of the traversal 
tree correspond to what we previously called forward connections. In the graphical 

29 



M'{t) =5t^ + 60t^2 + 1105 1^® + 27120 + 828250 + 30220800 
+ 1282031525 + 61999046400 + 3366961243750 
+ 202903221120000 + 13437880555850250 + 970217083619328000 1 
+ 75849500508999712500 + 6383483988812390400000 
+ 575440151532675686278125 + 55318762960656722780160000 
+ 5649301494178851 172304968750 
+ 610768380520654474629120000000 
+ 69692599846542054607811528918750 1^^^ 
+ 8370071726919812448859648819200000 + o(t^^°) 

Table 12. Development of the generating series M'{t), up to order one 
hundred and twenty. It gives as the coefficient of t^" the number of connected 
unrooted unlabeled triangular maps with n arcs, thus n/2 undirected edges 
and n/3 triangular faces (A062980). If we denote by a„ that coefficient, 
the recurrence is as follows: ai = 5 and for n > 1, a„+i = (6n + 6) a„ + 

En— 1 



Msit) = 3 1*^ + 11 1^^ + 81 1^^ + 1228 1^^ + 28174 1^^ + 843186 1^^ + 30551755 

+ 1291861997 1^^ + 62352938720 + 3381736322813 

+ 203604398647922 + 13475238697911184 + 972429507963453210 1 

+ 75993857157285258473 t^^ + 6393779463050776636807 

+ 576237114190853665462712 + 55385308766655472416299110 1^"^ 

+ 5655262782600929403228668176 

+ 611338595145132827847686253456 1^^^ 

+ 69750597724332100283681465962492 t^^o ^ ^(^^120^ 

Table 13. Order one hundred and twenty development of the generating 
series Ms^t) giving the number of connected unrooted unlabeled triangu- 
lar maps with n arcs, thus n/2 undirected edges and n/3 triangular faces 
(A129114). 



representation, we use an embedding of the produced triangulated polygon in the 
Poincare disc model of the hyperbolic plane as it seems the natural setting for generic 
non-overlapping triangular tilings. The surrounding circle around each figure is of 
course irrelevant to the structures. 



The tables 12 and 13 give the number of rooted triangular maps and unrooted 
triangular maps in the form of generating series. They are also part of [20j under 
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Table 14. The number of rooted triangular maps by genus (horizontally) 
and number of faces (vertically). 
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Table 15. The number of unrooted triangular maps by genus (horizontally) 
and number of faces (vertically). 



the references (A062980) and (A129114). Their computation is very similar to that 
of tables [6] and [7| which is explained in detail in [2l] . We shall explain in [23] in 
a unified fashion how one can compute generating series for rooted and unrooted 
unlabeled maps of various kind and in [22] the unexpected relation of this sequence 
to the asymptotcis of the Airy function. 

As another byproduct of the exhaustive list obtained from the generating algo- 
rithm, one can get the precise number of rooted and unrooted triangular map having 

summarize 



a given genus and a given number of triangular faces. Tables [14| and 15 
those results for small number of faces. Recently, M. Krikun [14] kindly communi- 
cated us recurrence relations satisfied by the entries of table [14] which he obtained 
by a clever recursive decomposition of rooted triangular maps. Those recurrence 
relations make it possible to evaluate easily those numbers without running the 
generator algorithm. Unfortunately, no general recurrence relation is known for the 



entries of table 15 The thesis [251 and articles [28l |29] contain the first enumerations 



of rooted maps of a given genus. 

Some lines of those two tables were previously known. For instance, the first line 



of table 14 is the number of spherical rooted triangular maps by the number of its 



faces [TS]. The first line of table 15 is its unrooted counterpart. It is computed by 
impressive closed formulae in a recent paper by Liskovets, Gao and Wormald [16] . 
The diagonal terms of those two tables also received close attention. For instance, in 
[To] Harer and Zagier computed the Euler-Wall characteristic of the mapping class 
group of once pointed genus g closed oriented surfaces by a remarkable combinatorial 
reduction of the problem in which rooted combinatorial maps with one vertex are 
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counted by genus yielding the diagonal sequence of the first table 1, 105, 50050, 
56581525, .... The diagonal sequence of the second table: 1, 9, 172, 1349005, .... 
gives the number of unrooted triangular maps of genus g with only one vertex. It 
has been studied at depth in the article [1] by A. Vdovina and R. Bacher. 



6. Concluding Remarks and Prespectives 

The generating algorithm presented in this paper (section [s]) may lead to trivial 
adaptations to generate wider classes of diagrams and combinatorial maps, possibly 
with prescribed degree lists for vertices or faces. Basically, it can be simply gener- 
alized to produce any connected pair of permutations with prescribed cyclic types, 
up to simultaneous conjugacy. 

Another way to extend the study, would be to modify the Dispatch procedure 
of section [s]) to generate not an exhaustive cover of the partial cases, 
but instead a single case of them picked at random. This would result, in a fairly 
straightforward fashion, in a random sampler algorithm of the corresponding com- 
binatorial structures instead of an exhaustive generator. The difficulty there, is to 
precompute precise conditional probability tables in order to control the probabil- 
ity distribution of the generated structures by bayesian techniques. Such tables of 
conditional probabilities could be computed with the help of generating series tech- 
niques namely by following the particular recursive structure of the algorithm and 
translating this recursive structure in functional equations on the generating series. 
This appeal for a further investigation and could be dealt with in a subsequent 
paper. 

If one designs the Dispatch procedure to pick at random with uniform 

probability this would unfortunately not result in a uniform distribution of the out- 
put structures but it can still prove useful to produce test cases to many algorithms 
operating on triangulations. A much beter approche but still not perfect is to gen- 
erate at random the two permutations a, and cTq having the right cycle types with 
uniform distribution and rejecting each time the resulting map is not connected. 
This procedure produces each map with probability proportional to 1/A where A is 
the number of its automorphisms. Since the maps have no automorphism other than 
identity in the vast majority of cases, that statistical bias is in practice unnoticeable. 
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